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C^) "Research on thermoelectrical energy conversion, the reuse of waste heat produced by some mechanical or chemical processes 
CN to generate electricity, has recently gained some momentum. The calculation of the electronic parameters entering the figure 
O of merit of this energy conversion, and therefore the discovery of efficient materials, is usually performed starting from the 
D .Landauer's approach to quantum transport coupled with the Onsager's linear response theory. As it is well known, that approach 
-suffers of certain serious drawbacks. Here, we discuss alternative dynamical methods that can go beyond the validity of the 
Xandauer's/Onsager's approach for electronic transport. They can be used to validate the predictions of the Landauer's/Onsager's 
approach and to investigate systems for which that approach has shown to be unsatisfactory. 
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;1 Introduction 

Electrical energy is one of the most versatile forms of energy, 
since it can be easily stored, moved, or converted to any other 
form (work, heat, etc.). On the other hand, in the process of 
electrical energy production, a relevant part of the energy, in- 
deed as much as 60% in most cases, is wasted as heat, e.g., 
the plumes from the cooling towers of a power station. The 
efficiency of this conversion process is limited by fundamental 
laws and, nonetheless, it is difficult to do much better than the 
present values. Nonetheless, no physical principle forbids us 
from using part of the waste heat of the first energy conversion 
to obtain more energy. Indeed, if the conversion process were 
ideally efficient no energy could be obtained through any other 
physical or chemical process. In reality, for non ideal cycles, 
by reusing the waste heat, we will pay a price in wasted energy, 
but looking at the overall process we will end up producing 
more electrical energy and less waste heat. This two-step pro- 
cess is limited in efficiency by thermodynamical laws, but it 
is still more efficient than the single step process. Many alter- 
native methods to convert that waste heat back to energy have 
been proposed. We will focus on one of them, i.e., thermoelec- 
tricity, the direct conversion of a temperature gradient into an 
electrical current.—— 

Nanoscale devices have the power to revolutionise the field 
of electronics. 6-9 It is clear that their comprehensive and ex- 
haustive experimental exploration not driven by reliable theo- 
retical predictions will be the result of serendipity not of scien- 
tific research. In the last few decades, Numerical and analytical 
tools have been developed to help understanding and designing 
efficient molecular devices. Remarkable successes have been 
scored, but the same tools for nearly all the most promising 
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molecules, e.g., oxides for photovoltaic devices, have given 
mixed results. The general consensus goes in the direction 
of investigating the same devices with dynamical methods.— 
These allow a careful description of the dynamics of both the 
electronic and atomic parts of the device thereby understand- 
ing their interplay in establishing the electrical properties of 
the molecule. One aspect that is attracting increasing interest 
is how to make these molecular devices more stable. We know 
that thermal excitations either destabilise the chemical bonds 
that attach the molecule to the leads or effectively destroy the 
lead by breaking the metal-metal (usually Au) chemical bonds 
of its components. 11 The investigation on how to efficiently re- 
move the heat locally stored in the device is thus of primary im- 
portance towards a nanoscale electronics.— Here again, time- 
dependent methods are the tools of choice, since they allow a 
fundamental description of the approach to a steady stated— 
(when a steady state is established) than the standard linear re- 
sponse theories. Finally, nanoscale devices seem to overcome 
some of the thermoelectric efficiency problems highlighted in 
bulk materials. A careful choice of the materials and an accu- 
rate tailoring of their properties may pave the way to efficient 
thermoelectric nanodevices. 

The efficiency of the thermoelectric energy conversion is 
related to the so called figure of merit. This number is defined 
in terms of the microscopical quantities of the devices, i.e., 



ZT^T 

K 



(1) 



where S is the Seebeck coefficient, a the electrical conduc- 
tance, k the thermal conductance, which includes the electronic 
and crystal contributions, and T the operating temperature. The 
figure of merit determines how efficient the Peltier or Seebeck 
devices will be: 
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where r) is the efficiency and r/o is the efficiency of the Carnot 
ideal cycle, r/o = (T> — T , < )/T > . Here, T> is the temperature 
of the hot junction, while T< is the temperature of the cool 
junction. Typical values for ZT in actual devices are around 1, 
and viable devices for technological applications need to have 
ZT larger than 2-3 M 

Nanomaterials have shown a plethora of interesting phe- 
nomena, many of which in the realm of electric and thermal 
transport 8 '2*i2.. For application to the thermoelectric technol- 
ogy, certain nanomaterials like metal nanowires present the in- 
teresting feature that electrical and thermal transport are al- 
most independent, in particular the electrons are mostly re- 
sponsible for the electrical transport while lattice vibrations 
are mostly responsible for thermal transport. This allows us 
to tailor nanoscale devices in which the thermal conductance is 
minimised while keeping the electrical conductance constant. 
On the other hand, quantum mechanical effects can strongly 
modify the transport properties of nanomaterials. As discussed 
by Hicks and Dresselhaus, 19 a key characteristic for determin- 
ing the efficiency of the thermoelectric energy conversion is 
the electronic density of states of the device. In small strongly 
confined systems, a tailored density of states might induce a 
manyfold increase of the figure of merit. For this reason there 
is widespread belief that nanoscopic devices can show a large 
thermoelectric figure of merit. 19 There are other reasons for 
searching for efficient thermoelectric materials in sub-micron 
world. We have seen that for large energy production thermo- 
electric devices might prove not efficient enough. However, 
their efficiency remains almost constant independent of the size 
of the device itself, while the efficiency of other energy conver- 
sion processes scales dramatically with size. For this reason, 
in applications in which controlling the temperature of a small 
portion of a larger device (like for example in cooling down a 
CPU) is fundamental, thermoelectric devices offer an economi- 
cal and efficient solution. In the quest for efficiency, the proper- 
ties of a wide range of materials have been investigated. Many 
have shown good thermoelectric efficiencies and this implies 
that in choosing the material for a certain application other pa- 
rameters can be adjusted, for example, it is easily foreseen that 
silicon nanowires might found wider implementation in current 
electronic devices based on silicon. Opaque or transparent ma- 
terials might find application in the field of photovoltaics, to 
increase the total current output of the device by using the part 
of the solar energy that is lost to waste heat. 

Many strategies have been proposed to increase the figure 
of merit of a thermoelectric device: we can either increase S, a, 
or the power factor S 2 a or decrease the thermal conductance k. 
It should become apparent that the maximisation of the figure 
of merit in bulk materials is difficult. Indeed, it is well known 
that the electrical and thermal conductances in bulk materials 
are not independent, but, for example in metals, they are con- 



nected by the Wiedemann-Franz empirical law 20 Q This law is 
based on the observation that in a metal, electrical and thermal 
currents are both carried by electrons. However, little is known 
about how to calculate the Seebeck coefficient. The definition 
of S is quite generic, 



where AV is the voltage drop at the junctions if a temperature 
gradient AT is maintained, in an open circuit configuration, 
i.e., no electrical current is flowing. It has been suggested that 
a direct measure of this coefficient gives direct evidence of the 
sign of the carriers, a problem in complex nanoscale systems.— 
This statement is generally true for bulk materials, but it has 
been pointed out that the theoretical approach used to under- 
stand the experimental results is based on the assumption that, 
at equilibrium, both the thermal and the electrical currents van- 
ish. 22 This assumption is not valid when we consider nanoscale 
systems, i.e., the ones for which we hope to get large figures of 
merit. It has been shown that removing this condition allows 
different signs of the Seebeck coefficient given the kind of car- 
riers in the device while varying the working conditions.— 

In this Perspective, I will discuss a nowadays standard ap- 
proach to the calculation of the figure of merit, namely the com- 
bination of the Landauer's approach to quantum transport and 
the Onsager's theory of non-equilibrium response. I will pro- 
pose alternative dynamical methods that can go beyond certain 
limitation of that standard approach. The paper is organised 
as follows. In the following section, we will discuss the Lan- 
dauer's theory of quantum transport. We will show how one can 
arrive to simple expressions for the physical parameters enter- 
ing the figure of merit ZT. In section [3] we will discuss an 
alternative method for the calculation of the electrical conduc- 
tance. In section |4] we will turn our attention to the Seebeck 
coefficient. Finally, in section we will discuss the problem 
of evaluating the thermal conductance for electrons from a dy- 
namical theory. In this latter case, we cannot provide an answer 
to this question, but we will discuss briefly two research lines 
which could lead to a comprehensive dynamical theory of the 
thermoelectric energy conversion. Throughout the paper we 
will set H = &b = 1, where ks is the Boltzmann constant. 

2 Landauer's approach to quantum transport 

A widely used and successful approach to the calculation of 
the conductance in nano- and mesoscopic systems is the Lan- 
dauer's formula. 18 In this elegant and economic formalism one 
assumes that the details of the system are irrelevant far from 



* Notice that the Wiedemann-Franz law could be formulated in terms of conduc- 
tivities or conductances. For our scope, this distinction is irrelevant, since the 
same dimensional factors will multiply the numerator and denominator of[T] 
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the nanostructure in some region of space called reservoirs: the 
only property requested from these reservoirs is that their cross 
section be infinite, meaning that an electron entering that region 
cannot escape. A fixed bias keeps the two reservoirs polar- 
ized. The nanostructure is connected to the reservoirs through 
two regions of space which smoothly enter it. The transport 
properties are completely determined by the probability am- 
plitude that an electron comes out of a reservoir, crosses the 
lead-nanostructure-lead device and is absorbed into the other 
reservoir. It is clear that this model contains many assump- 
tions. We should be able to identify the different regions of in- 
terest (reservoir-lead-nanostructure) where different physics is 
at play, a separation lacking a solid justification from a micro- 
scopical point of view.— Another basic assumption is that the 
many-particle states in the reservoirs are fully determined and 
do not mix; we are effectively neglecting particle-particle in- 
teraction. Finally the Landauer's approach cannot explain the 
nature of out-of-equilibrium states that are established in the 
device. In particular, it cannot answer the basic question, will 
such a system ever reach a steady state? Other approaches suf- 
fer from similar drawbacks. For example the non-equilibrium 
Green's function formalism (also known as the Keldysh formal- 
ism), 24 although a powerful technique, assumes the existence 
of a steady state imposed by scattering boundary conditions far 
from the nano-system. The existence of a steady state seems 
a reasonable assumption for mesoscopic systems but needs to 
be proven for nanoscopic devices or closed system. 16 Other ad- 
hoc models have been proposed in the past but they lack a clear 
physical interpretation.— 

Assuming the Landauer's approach can be applied, we start 
from the linear out-of-equilibrium Onsager's relations for the 
currents flowing along the deviceQ 



j = L 00 Afi + L iAT 
j h = L w Afi + in AT. 



(4) 



Symmetry considerations lead to Lqi = -^io- Landauer's ap- 
proach to thermal transport allows to express the Onsager's co- 
efficients L's in terms of the transmission probability. In linear 
response with respect to the thermal gradient AT and the po- 
larisation potential A/i, one obtains the general expressions 



Lij{p.,T) 



/CO 
dE(E - 
-CO 



df(E,T,n) 
dE 



T{E). (5) 



In (0, the chemical potential /i = (/ir + /il)/2, where [1r and 
Hl are the chemical potentials of the reservoirs, respectively. 
By assumption, — Hr = A/x <C /ir. The same is true 
for the temperature T = (T R + T L )/2, with AT = T L - T R < 
Tr, Tl. f(E,T,fi) is the equilibrium Fermi distribution. We 

f Here the particle and heat currents are assumed to be perpendicular to the cross 
section of the device at each position. We expect this to be correct for long 
nano-wires. 



assume indeed that the reservoirs are large and at equilibrium 
at given temperature and chemical potential. Moreover it is as- 
sumed that the electrons in the reservoir do not interact. This 
last assumption is justified by the fact that in metal the screen- 
ing length is small. Finally, T(E) is the transmission proba- 
bility for one electron with given energy E to travel across the 
device. The theory needs to calculate this quantity to make the 
set of equations (0]i and (0 of any use. We would like to point 
out here that while Landauer's approach to transport is valid 
beyond linear response, the Onsager's relations are valid only 
in linear response. 

Given the Onsager's relations, and remembering the defini- 
tion of the different physical quantities, we obtain the follow- 
ing relations between the microscopic parameters that enter the 
figure-of-merit and the integral L's, 



S 



e L Q0 (fj,,T), 
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In (0, K e is the electronic contribution to the thermal conduc 
tance, defined by 

rri Jh 

K e l = — — - 

AT 

^ 3=0 

Here, we assume that we can write the total thermal conduc- 
tance as k — n e + Ki, where ki is the crystal (either phononic 
or vibronic) contribution to the thermal conductance. This is 
valid if the electron-phonon interaction is small. Only then we 
can separate these two contributions. 

A step further into a complete theory of the thermoelectric 
energy conversion is the expression of the thermal conductance 
ki in terms of the transmission probability of a phonon across 
the device, 



Kl = 



dT 



(8) 



where n(ui,T) is the Bose distribution of non-interacting phonon 
at given temperature T, and t(cj) is the probability that one 
phonon with energy u> could cross the device.— ~— 

It is clear from (0, ((6), and ^ that the only ingredients that 
we are still missing are the transmission probabilities T(E) and 
t(lo). A lot of progress has been made to connect these prob- 
abilities with the non-equilibrium Green's function theory and 
ultimately with electronic structure calculations. 27-30 Remark- 
ably, the transmission probability T can be written as 



T(E) = TT[G r (E)T L G a (E)T R ] 



(9) 



where G r is the retarded Green's function, G a (E) = G r {E)^ 



and T L 
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J L(R) 
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describe the coupling, via 
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the self-energies T, r and S a of the left and right reservoir with 
the device, respectively. Finally, the non-equilibrium Green's 
function is obtained from the Hamiltonian via, 

GT{E)= [(E + iQ+)+H c -^ r L -W R \- 1 , (10) 

where He is the Hamiltonian of the central system. The On- 
sager's out-equilibrium response theory is exact in the linear 
approximation. Therefore, the quality of the results is strictly 
related to the approximations used to calculate the Hamiltonian 
He and its eigenvalues. A certain degree of success has been 
obtained by using for the central region the Kohn-Sham (KS) 
Hamiltonian. 18 i 31 i 32 To some extent, this degree of success is 
unexpected. The KS theory indeed, should not give access to 
the spectral properties of the system under investigation. Mean- 
while, those properties are indeed strongly used in the evalua- 
tion of the non-equilibrium Green's functions. Finally, expres- 
sions similar to (O and < TT~0b can be written for the phonon- 
transmission probability, with the substitution of E with w 2 , 
He with K the spring constant matrix, G r with D r the phonon 
non-equilibrium Green's function. 

The fact that one can tackle on the same ground both the 
vibrational and electronic contribution to the figure-of-merit is 
one of the greatest advantages of this theory. Moreover, its sim- 
plicity in expressing all the quantities of interest in terms of the 
transmission probability and this one in term of a microscopic 
model of the device, is invaluable in making the theory predic- 
tive.^- 

At the same time, to gain a deeper understanding on the 
physics of thermoelectric energy conversion, we should be able 
to build other models and theories that can extend or confirm 
the validity of the Landauer/Onsager's approach to the calcu- 
lation of the figure-of-merit. Over the years, many novel ap- 
proaches have been proposed which are trying to go beyond 
the actual state-of-the-art in molecular transport, by, e.g., im- 
plementing a GW approximation for the calculation of the non- 
equilibrium Green's function and apply it to transport in molec- 
ular device s 34,35 or by building a suitable base for calculating 
the time evolution of the electrons. 36 In the following we will 
discuss other ways, strongly based on Time Dependent Density 
Functional Theory (TDDFT), 10 to calculate the Seebeck coef- 
ficient and the electrical conductance by using dynamical the- 
ories. Here, one studies the dynamics of the electrons as they 
respond to the external perturbation that generates the electrical 
and heat flows. In general, these theories have the possibility 
to go beyond the linear response, and assess non-linear effects. 
The drawback is that the computation of the electron dynamics 
is more expensive than any static calculation. At the same time, 
the dynamical theories are still in their infancy. This means 
that we do not have reliable approximations to the many-body 
Hamiltonian that allow a correct evaluation of the dynamics of 
the observable, although some remarkable progress has been 
made through, for example, exact solutions— or in correlating 



the dynamics of the electrons with the ions. 38,39 In this respect, 
further investigation is necessary. While these alternative mod- 
els do allow the evaluation of the Seebeck coefficient and the 
electrical conductance, little can be said about the thermal con- 
ductances, especially for the electrons. In Section [5] we will 
shortly address this question and outline two possible research 
paths for modelling the thermal electronic conductance with 
dynamical theories. 

3 Time-dependent (current-)density functional 
theory 

Density Functional Theory (DFT) has become a ubiquitous tool 
in the investigation of the properties of matter.^ 3 - 2 - The reasons 
for this success are manifold and introducing the reader to the 
field is well outside the scope of this article. It will suffice to 
say that DFT gives in principle access to an exact solution of 
the many-body problem, whenever we can focus only on cer- 
tain quantities, like, e.g., the single particle density. In its origi- 
nal formulation, DFT deals with the problem of calculating the 
ground state energy and single particle density of a many-body 
electronic system. Over time, many extensions of the original 
idea have been presented: the most interesting for our discus- 
sion are the TDDFT and the Time-Dependent Current-Density 
Functional Theory (TDCDFT). 10 i 32 i 40 i 41 As their names sug- 
gest, the two theories deal with the problem of calculating the 
dynamics of a many-body system. The result is obtained by 
mapping the dynamics of the real system onto the dynamics 
of a fictitious many-body system of non-interacting particles. 
In the case of TDDFT, we are interested in the single-particle 
density, n(r, t). Interestingly, TDDFT does not give access to 
the current density j(r, t)— but it can reproduce the exact total 
current.— On the other hand, TDCDFT gives access to the ex- 
act single-particle current density, and via the continuity equa- 
tion, to the single-particle density. In this respect TDCDFT is 
a more general theory than TDDFT. We would like to point out 
that DFT and TDDFT should be able to deal with the dynamics 
of strongly correlated systems. In this direction, remarkable re- 
sults have been obtained in the solution of the Kondo problem, 
see for example 44 and references therein. One of the advan- 
tages of this theory over the standard non-equilibrium Green's 
function formalism is that it avoids the separation of the sys- 
tem intro three (or more) parts. In general, indeed for the non- 
equilibrium Green's function we need to identify the left and 
right reservoirs and connect them with the central region. In 
most of the cases, this procedure lacks a clear physical inter- 
pretation 23 and may produce spurious effects. On the other 
hand, in TDCDFT one can investigate the full dynamics of the 
system without reverting necessarily to this separation. 13 i 45 i 46 

It should also be clear that TDCDFT is built to treat prob- 
lems where electrical transport is important. By reproducing 
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coupling has the form 

H coup = \V®B, (13) 

where V and B are system and environment operators, respec- 
tively, and A is a small coupling parameter, a standard proce- 
dure leads to the equation of motion for the "state'H | "J) of the 

syste m 52 i 54 ^ 

i^ t mt))=-iH C \m) 

+ A 2 f dTC(T)V^T)e- iHcT V\^(T-t)) 
Jq 

+ \r,(t)V\*(t)), 

(14) 

where C(t) = Tr E [pf B(t)B^(0)] is the environment corre- 
lation function. We assume that the environment is maintained 
in the thermal equilibrium described by the statistical operator 

P " = Tr E e-P"' ' (15) 

and rj(t) is a complex coloured noise with the statistical prop- 
erties 52 



the single-particle current for given external potential, the elec- 
trical conductance, a can be obtained directly from its micro- 
scopical definition 

j(r,t) = J dr'a(r,r') ■ E(r',t). (11) 

Since the electrical field E(r,t), is an external input, (TTTt al- 
lows the complete determination of the conductance tensor, if 
the current density can be calculated. Moreover, the theory 
does not suffer from any problem of dividing the device/reservoir 
region. In principle, we can treat the motion of the electrons ex- 
actly in the presence of an external polarizing potential.— 

4 Calculation of the Seebeck coefficient 

We will now turn our attention to the calculation of the Seebeck 
coefficient. It should be apparent from its definition (j3]l that the 
calculation and measurement of the Seebeck coefficient are a 
difficult task. Indeed, in the experiments we combine the abil- 
ity to control the temperature gradient with the measurement of 
small bias drops while maintaining an open circuit configura- 
tion. These technical difficulties reflect into a large uncertainty 
in the values of the Seebeck coefficient which is usually the re- 
sult of indirect measurements which in turn require a good level 
of theoretical understanding. From a theoretical point of view 
is then of the uttermost interest to have solid and reliable ways 
to calculate the Seebeck coefficient for many different materi- 
als. In this respect, the theory we have exposed in Section |2] 
has proven an invaluable tool. Nonetheless, we must find an al- 
ternative formulation that will teach us the strengths and weak- 
nesses of the Landauer's approach and provide novel tools to 
pursue our quest of efficient thermoelectric materials. Ideally, 
we would like to revert to the basic definition of the Seebeck co- 
efficient, S = —AV/AT at vanishing electrical current. AV 
is the response of the system to the external thermal gradient, 
AT. We then need a theory able to describe the dynamics and 
steady state of a quantum mechanical system coupled to two 
external baths, or, more in general, in contact with an external 
environment. 

The theory of open quantum systems has been devised with 
this aim in mind.— ~— The general starting point is the Hamil- 
tonian of the composite system: device plus external environ- 
ment. Generally we can write that Hamiltonian as 

H = He + He + H coup (12) 

where He is the Hamiltonian of the isolated system, He the 
Hamiltonian of the environment, and H coup describes the cou- 
pling between the two. The aim of the theory is to "fold" the 
degrees of freedom of the environment and obtain an effective 
equation of motion for the system dynamics. If we assume the 



V (t) = 0, V (t) v (0) = 0, r?(t)r?*(0) = C(t). (16) 

Here with the 777 we indicate the average taken after many re- 
alisations of the noise. 

The stochastic Schrodinger equation (SSE) ( fT4l is able to 
reproduce the dynamics of the density matrix for the system^iS 



p=LJ==d. (17) 

(*l*> 

Interestingly, it is the average over the many realisations of the 
noise that makes this density matrix not pure. 

A standard approximation consists in replacing the correla- 
tion function C(t) with a Dirac S(t) function. This corresponds 
to neglecting the time correlation of the bath, effectively assum- 
ing the environment has no memory - a Markov approximation. 
Within the Markov approximation, the SSE assumes the simpli- 
fied form 

i^-Wt)) =-iH \9(t)) 
gt \ \ II oi v ii (lg) 

+ A 2 VV|*(t)> + Xr](t)V\^(t)}. 

Using the definition, we can derive a simple form for the dy- 
namics of the master equation, 

^ = -i [H, p] - VpV^ + \v^V P + \pV f V. (19) 

| It is important to point out that the state \t generally does not describe a true 
quantum trajectory^! 
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It can be proven that this is the most general memoryless mas- 
ter equation valid up to second order in A that preserves the 
positivity of the density matrix.— 

The equations (fl~8T > and ( [T9| > can be easily generalised to the 
case of the system being in contact with more than one bath 
operator B. In doing so, we can assume that each of these 
baths corresponds to a different thermal equilibrium (for ex- 
ample with different temperatures) and then study the dynam- 
ics of the system under this now out-of-equilibrium conditions. 
This idea has been explored partially to investigate the thermo- 
electric energy conversion with (fT9ll. 22 It is important to point 
out that in principle, in the Markov approximation, the detailed 
balance equation that will ensure the establishment of thermal 
equilibrium when the temperatures of the baths are the same, 
is not fulfilled. 54 Moreover, it is not clear how to treat particle- 
particle interaction in this formalism, since that in general will 
make the Hamiltonian state-dependent. The derivation of the 
master equation in this case appears more subtle and care is 
needed to avoid spurious effects. 60 For this reason, a TDCDFT 
scheme for this open quantum system has been recently devel- 
oped that can go beyond the limits of the density matrix for- 
malism.^ 6 ^ This scheme, by coupling the strength of a DFT 
approach with the theory of open quantum systems, could help 
in establishing a framework for the investigation of the See- 
beck coefficient of materials. Indeed, being able to couple the 
device of interest with two baths at different temperatures, we 
can study the formation of the potential difference AV in time, 
by using the Poisson equation, since the electrons in the system 
will try to oppose to the thermal gradient by piling up closer 
to the coolest bath. At equilibrium, no electronic macroscopic 
current will be present in the device, and we can use the defini- 
tion (0 to evaluate the Seebeck coefficient. 22 The same theory 
can be used to investigate the validity of the Fourier law at the 
nanoscale.— 

5 Thermal conductance of electrons 

As we have seen in the evaluation of the figure of merit, the 
Seebeck coefficient and the electrical conductance can be cal- 
culated by using dynamical methods thus offering an alterna- 
tive description to the Landauer formalism. 

The situation is less clear with respect to the calculation of 
the thermal conductance. In this case, there is not a dynami- 
cal formulation of the problem in terms of a density functional 
approach. The reason for this situation is simple: The energy 
current jh cannot be written in terms of the electrical (or den- 
sity) current. Indeed, we can define the energy current across a 
given section of the system as 

Jh = ~J t {Hs) (20) 
where Hs is the energy stored in the volume enclosed by the 



surface we are considering. For example, if we can separate 
the system into left, right and central region, where the left and 
right act as energy "reservoirs", Hs is the energy stored in one 
of these two reservoirs. We could then write a similar equation 
for the energy coming from the other reservoir. Combining the 
two, we can calculate the energy flow in the system. A naive 
approach would simply suggest to calculate the energy starting 
from the DFT wavefunctions, in a fashion similar to what has 
been done in the past to evaluate the non-equilibrium Green's 
function. However, this operation contains a certain number of 
approximations, many of which are not under control. Indeed, 
the theorem of DFT (either the static or dynamical) cannot give 
access to the energy of the excited states. For this reason, the 
concept of the total energy stored in a part of the system will 
bear little meaning in a DFT formulation of the problem. 

We foresee two ways out of this impasse. On the one hand, 
one can think of building an ad-hoc density functional theory 
that provides the exact energy or thermal current. To do so, 
we need to identify a suitable potential that is connected to 
the energy current, responsible to generate that current in the 
many-body system. If that potential is found, its existence is 
doubtful to say the least, then we can hopefully prove that the 
mapping between the potential and the energy current is one- 
to-one. If this can be done, and a suitable KS scheme be built, 
we will have an exact formulation of the problem in terms of 
non-interacting particles and we can compute the energy or the 
thermal current from first principles. 

However, one can build some approximate model to the 
electron thermal transport. In this respect a hydrodynamical 
formulation of the electronic transport, would provide a set of 
equations which includes the thermal conductance and the spe- 
cific heat of the electron liquid. See for example Ref.— for 
the formulation of the theory for a quantum point contact. In 
Ref. 12 the hydrodynamical equations were used to calculate the 
electronic temperature in the nanojunction. However, if we do 
assume that temperature can be either modelled or measured in- 
dependently, these equations will be sufficient to calculate the 
thermal conductance of the electrons. An ab-initio theory here 
will enter in defining the viscous coefficient of the electron liq- 
uid. 

In a hydrodynamic approach to quantum transport, we de- 
fine the stress tensor starting from the velocity of the fluid, 

v(r,t) = j(r,t)/n(r,t), 

7Ti,j = r) \ ViVj + VjUi - ^Sij'VkVkj (21) 

where 77 is a real coefficient (the viscosity) that is a functional of 
the density. 12 i 63 - 64 We point out that ( T2"T1 i is in fact a particular 
case of a general stress tensor with memory effects taken into 
account. 65-67 The parameters of the liquid can be build from 
a microscopic theory of the low-energy excitations of the sys- 
tem under investigation. Unfortunately, it turns out that their 
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detailed dependence on the geometry and dimensionality of the 
system is complex. 64 It is clear that further investigation in this 
direction is needed. In terms of this stress tensor, the continuity 
equation and the equation of motion for the velocity field are 
written as, 

D t n{r, t) = — ri(r, i)V • v(r, t), 
mn(r,t)D t Vi(r,t) = -ViP(r,t) + V j ir i j(r,t) (22) 
-n{r,t)ViV ext {r,t). 

where D t — d/dt+v(r, t)-V is the convective derivative, rn the 
electron mass, P(r, t) is the pressure, and V ex t (r, t) is a time- 
dependent external potential. The second equation is nothing 
else than the Navier-Stokes equation for the a viscous fluid.i2i£& 
It is interesting to point out that for a nanoscale system, for 
electronic low-energy excitations, this theory is equivalent to 
the many-body formulation. 

To study the thermal transport, we need to complement (1221 
with an equation for the energy flow. For a quantum point con- 
tact, a nanoscale constriction between two large metallic reser- 
voirs, one obtains,— 

Wi d (r)djViC) + v ' [fc(r)VT e (r)] = c v {T e )v{r) ■ VT e (r), 

(23) 

where T e is the electronic temperature, k(r) is the diffusion 
constant and cy is the specific heat at fixed volume of the elec- 
tron gas. Obviously, in writing Eq. d23l we have assumed that 
some thermodynamic quantities like temperature and entropy 
for an electron liquid flowing in a nanostructure can be defined. 
This is a much debated point, and obviously we do not have a 
general solution for it. However, here we argue that the elec- 
tron temperature may be defined as the one ideally measured 
by a probe weakly coupled to the system and in local equilib- 
rium with the latter.i&22 The set of hydrodynamical equations 
has been used to study the thermal heating of a nanoscale de- 
vice. They do predict novel and interesting phenomena, like 
a cooling of the phonon modes due to the electron heating.— 
This result has been beautifully confirmed by experience. 69 

I am personally amazed that the equation of motion for a 
many-body system does reduce to a form similar to the Navier- 
Stokes equations of a classical liquid. The quantum mechanical 
effects are still embodied by the pressure term P(r). For this 
reason, the solution of these Navier-Stokes equations is some- 
how more difficult than the already demanding solution of the 
classical equations. Finally, it is quite surprising that TDCDFT 
arrived at this same conclusion well before.^£& At the same 
time, TDCDFT seems to suggest natural mechanisms for en- 
ergy relaxation, 70 and recently a thermodynamical interpreta- 
tion of the KS energy has been given, as the maximal work that 
can be extracted from the system. 71 Starting from these consid- 
erations, it does appear possible that TDCDFT can be used to 
describe some of the physical process related to energy dissi- 
pation and transfer. 



Both formulations of the problem of calculating the elec- 
tron thermal transport will require the solution of long-standing 
problems. For example, in a density functional theory for the 
thermal current, we will very soon discover that a thermal cur- 
rent is generated by almost any external potential, not only by 
a temperature imbalance. Therefore, how can a given thermal 
current be in a one-to-one mapping with a certain potential? 
The hydrodynamical formulation can be equally as challeng- 
ing. For example, what is the range of validity of the equations 
of motion. For sure, at short wavelengths, the electrons do not 
behave as a liquid. Moreover, at the nanoscale, a few very fun- 
damental equations of classical heat transport, like the Fourier 
equation, are not valid. The analysis of thermal transport at the 
nanoscale thus should include a wider analysis on the validity 
of assumptions made on the dynamics of the system. 

We have left out the problem of phonon/vibron thermal 
transport, and the coupling between the electrons and the vi- 
brons. Presently, besides doing a full dynamics of the cou- 
pled electronic and atomic motions, we do not see an alterna- 
tive to the static calculations performed starting from (JHJ, and 
improving the calculation of the Green's function via pertur- 
bation theory. A lot of effort has been put into the calcula- 
tion of the phonon thermal conductance, for example in sili- 
con based nano-structures,- 7 -^ via either molecular dynamics 
or the Boltzmann transport equation. While these techniques 
have been proven useful in calculating the thermal conduc- 
tance, it also remains difficult to obtain full convergence of 
the dynamics especially when the system is maintained out of 
equilibrium. On the other hand, it is true that © is limited to 
the harmonic approximation for the phonon spectrum. This ap- 
proximation appears to be valid only for low energy vibrations- 
since the interaction parameter vanishes with vanishing phonon 
energy-which are also responsible for the largest contribution 
to thermal transport. 

6 Conclusions 

The investigation of the thermoelectric properties of materials 
poses interesting challenges because it does require the treat- 
ment, on equal footing, of electrical and thermal transport. So 
far the only theory that has succeeded in doing so is the lin- 
ear response theory developed by using the Onsager's relations 
coupled with the Landauer's approach to thermal transport. While 
successful for many systems, it is necessary to have alternative 
theories able to describe the same systems and go beyond the 
non-equilibrium linear response regime. In this article, I have 
discussed the state-of-the-art of one of the alternatives, namely 
a dynamical approach to the calculation of the thermoelectric 
coefficients. While the theory is based on solid ground for the 
Seebeck coefficient and the electrical conductance, the ther- 
mal conductance, both for the electronic and atomic systems, 
is more difficult to access with first principle theories. I have 
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